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AGCELERATFi)  LIFE  TESTING  - 
A PR/\GM/\TrC  BAYESIAN  APPROACH 

by 
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0 . Introduction  and  Overview 

In  practice,  long  life  items  are  often  subjected  to  larger 
than  normal  stresses  (doses)  in  order  to  obtain  failure  data  in  a 
short  amount  of  test  time.  It  is  common  to  test  the  items  at  more 
than  one  stress  level,  and  typically,  more  items  are  tested  at  tlie 
higher  stress  levels  tlian  those  at  the  lower  stress  levels;  such 
tests  are  known  as  ac^ie legated  or  overs tresr,  life  tests.  The  basic 
aim  is  to  make  Inferences  about  the  life  distribution  of  the  items 
at  the  normal  stress  levels  using  the  failure  data  from  accelerateii 
tests. 

The  current  approach  to  this  important  problem  makes  inferences 
under  parametric  assumptions.  This  may  be  valid  in  some  situations, 
and  does  yield  results  which  are  appealing  from  a statistical  point 
of  view.  However,  an  engineer  or  a statistician  working  under  less 
well-defined  conditions  may  find  this  approach  too  restrictive.  For 
example,  it  is  common  to  assume  that  at  all  stress  levels,  failure 
times  are  governed  by  exponential  distributions  or  Welbull  distribu- 
tions. In  addition,  a functional  relationship  between  the  parameters 
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of  the  lailure  distribution  and  the  applied  stress  is  assumed.  Such 
a relationship  is  known  as  an  aijccUu'at  ion  function  or  a time  tconii- 
jotvution  fu>ictton;  examples  can  be  found  in  Mann,  Schafer,  and  Sing- 
purwalla  (Chapter  9,  (1974)).  In  two  recent  papers,  Phatak,  Zinuner, 
and  Williams  (1977)  and  Shaked,  Zimmer  and  Ball  (1977),  the  distribu- 
tional assumption  is  dropped,  but  the  requirement  is  retained  that 
tl»e  (unknown)  failure  distribution  be  of  the  same  form  at  the  use  and 
at  ail  the  accelerated  stress  levels. 

Clearly,  tliere  are  many  situations  where  the  above  assumptions 
may  not  be  appropriate.  Of  particular  concern  is  the  assumption  iliat 
the  failure  distribution  is  of  the  same  form  at  ail  the  stress  levels. 
One  reasi)n  for  this  concern  is  that  different  stress  levels  may  have 
different  effects  on  the  mechanism  which  causes  failure,  and  thus  from 
a physic. tl  point  of  view,  it  may  bo  more  realistic  to  allow  for  differ- 
ent forms  of  the  failure  distribution  at  tlie  different  stress  levels. 

The  approacti  that  we  propose  in  this  paper  requires  neither 
dlstributlon.il  assumptions  nor  the  specification  of  a time  transform.i- 
t ion  function.  Rather,  our  .approach  is  Bayesian,  and  is  prompted  by 
what  is  .ictually  done  in  practice.  The  Bayesian  point  of  view  .allows 
us  to  Incorpor.ile  some  j'n.'o!':'  Information  which  is  available  in 

accelerated  life  tests.  We  would  like  this  paper  to  be  construed  .is 
prellmln.iry  .in.l  pr.iy.matic,  and  thus  have  not  attempted  to  give  a t>ill 
theoretical  justification.  Consequently,  we  would  like  to  invite 
mathematical  statisticians  to  resolve  the  statistical  problems  posed. 

1.  Trc^l  imi_n;^r  ies 

As  stated  hetore,  in  accelerated  life  tests  the  items  are  tested 
under  different  stress  environments.  A stress  environment  may  be 
characterized  by  a single  stress  such  ns  voltage  or  temperature,  or 
by  multiple  stresses,  each  of  a different  type.  We  denote  a stress 
environmt'nt  by  F.  . and  the  sot  of  all  F's  by  ^ . We  assume 
that  the  elements  ot  may  be  ordered  according  to  the  magnitude 
of  their  severity.  Thus,  for  .iny  two  elements  Kj  and 
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belonging  to  ^ ^ denotes  the  fact  that  E^  is 


than  E 


j ■ 


Ls  more  severe 

Let  us  denote  the  k accelerated  stress  environments  by 


, E2 , . • . , E|^  , and  the  normal  or  use  conditions  stress  environment 


by  E 


We  shall  assume  that  the  set 


r i 


is  completely  ordered  with 


respect  to  the  relationship  ^ and  that 


E,  > E„  > . . . > E,  > E . 

12  k u 

The  basic  problem  in  accelerated  life  testing  is  to  make 
inferences  about  the  life  distribution  of  the  item  under  use  environ- 
ment E^  using  failure  data  obtained  under  accelerated  stress  environ- 
ments E^,  E^ , . . . , E|^  , where  k 1 • In  some  situations,  it  is  possible 

to  obtain  a limited  amount  of  failure  data  under  the  use  conditions 
environment  E^  ; however,  in  practice,  these  situations  do  not  appear 
to  be  very  common.  Although  having  failure  data  under  the  use  condi- 
tions environment  E^  has  advantages,  its  absence  is  in  no  way  detri- 
mental to  our  approach.  However,  in  order  to  obtain  results  wliich  are 
useful,  it  is  necessary  tliat  k be  moderate  to  large,  and  this  is  what 
we  require.  This  requirement  does  not  Impose  any  practical  difficul- 
ties in  many  situations  of  interest,  especially  those  involving  the 
accelerated  life  testing  of  electronic  components  and  in  bioassay 
experiments  on  animal  populations. 


In  order  to  introduce  some  notation,  let  us  denote  the  failure 
distribution  of  the  items  wliich  are  tested  under  environment  E^  by 

Fj  , where  (0  ) = 0 for  all  values  of  j . We  assume  that 

is  absolutely  continuous  and  thus  f^  (x)  , its  probability  density 

function,  exists  for  xg[0,“>)  . If  we  denote  1 - I’V  Ej(x)  » 

then  Xj  (x)  , the  failure  rate  of  Fj(x)  , is  defined  by 


f.Cx) 

X.(x)  = J , F (x)  > 0 . 

J F (x)  ■' 


3 


F 


T-  .57  . 


We  find  it  convenltMU.  .ind  reasonable  to  assume  that  V^(x)  is 

continuous  in  x for  xc(i),®)  . We  use  the  ti'rms  fit.' 

and  /ujaci.'H!  interchangeably.  .Since  Kj  >•  ^ •••  ^ ^ • 

it  is  logical  to  assume  that 

\ (x)  > X^fx!  •••  ^ X (x)  ^ X (x)  (1.1) 

for  all  xelO,-")  . 

Using  failure  data  obtaimxl  under  , . . . , 1’^^  , we  wouKl 

like  to  obta  in  estimators  (^)  ” l,2,..,k  , such  that 

for  some  0^1.*^"  and  all  xe{ 0,1.1  • 

8t  , St  St  , 

Xj^(x)  ^ X^Cx)  > •••  > X|.(x)  . n .2) 

St 

The  notation  .\  ^ Y denotes  the  fact  that  X is  stochast  ic;.  1 ly 

larger  than  Y ; i.e., 

P[X^xl  ^ P[Y^x]  lor  all  values  of  x . 

In  order  to  obtain  estimators  of  ’ -i  “ l,2,...,k  , 

which  satisfy  Equation  (1.2),  we  shall  use  a Bayesian  approach.  Under 
this  approach.  Condition  (1.1)  is  Incorporatei  as  a prior  assumption. 

Our  approach  is  in  contrast  with  that  of  Brunk,  Franck,  Hanson  .»nd 
Hogg  (1966),  who  embody  a similar  but  weaker  conditioit  than  (l.l)  in 
their  likelihood  function.  Specifically,  Brunk  et  al.  assume  that 
the  distributions  are  ordered,  i.e., 

F (x)  > F (x)  > •••  > F (x)  Vx  . 

1 — 2 — — k 

The  above  condition  is  a consequence  of  Condition  (1.1),  but  is  not 
equivalent  to  it.  Also,  the  approach  of  Brunk  et  al . is  not  Bayesian. 

We  shall  first  present  a methodology  for  a Bayesian  estimation 
of  the  individual  hazard  rate  Xj(x)  unconstrained  by  (1.1). 
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2 . Bayesian  Estimation  of  a Single 
Hazard  Rate  Function 


For  estimating  an  individual  hazard  rate  (t)  , we  use 
the  following  Bayesian  procedure. 

Let  number  of  items  exposed  to  the  environment 

K,  at  time  t . Thus,  N (0)  Is  the  number  of  items  tliat  are  initially 
j J 

put  on  test  in  environment  E.  . 

J 

We  hypothesize  that  the  failure  times  are  governed  by  a time- 
dependent  Poisson  process  with  the  probability  of  failure  in  [L,t+h] 


given  by 


Nj (t)Xj (t)h  + o(h)  , 


where  X (t)  is  the  failure  rate  at  time  t under  environment  E . 

J 

For  purposes  of  analysis,  we  divide  the  time  interval  [o,L] 
into  intervals  of  length  h > o , where  h is  chosen  to  make  L a 
multiple  of  h . For  convenience,  we  denote  (i-l)h  by  t^  , 

N (t  ) by  N , and  X (t  ) by  X , for  i = 1 , 2, . . . , (L/h)  . 

J i d.i  j i j.i 

Let  X.  , denote  the  number  of  failures  in  [t,,t.+h)  : if 
J ,1  i 1 

there  are  no  withdrawals,  removals,  censoring,  etc.,  then 

N,  , , * N,  . - X.  , . It  is  helpful  to  clarify  the  above  terms  by 

j.i+1  j.i 

the  following  diagram: 


Number  surviving  = N, 


lh=t2  2h=t2 


Number  of  failures  = x.  . 

.1.1 


(i-l)h»'t^  lh=t^^j^ 


1^1 


Ttie  loiiU  li  ist  r ihut  ii)n  ot  x.  ,,x,  ,....,x. 

i.l  j.J  + l 

given  by  the  multinomial  distribution 

(l./h)  + l X j 

'■/"jo'  • 


(2.1) 


where  p.  is  the  probability  of  failure  of  a specified  unit  in 
J * ^ 


(ti,ti^l)  . with 

(I./h) 

Pj.(L/h)+l  - 1 - • P.  , . ••‘nd  X 


i-1 


i.i 


(l./h) 

),(I,/h)  + l ’'j.i 


i-1 


Donot ing 


11  - - Pj  1 ’ (2.1)  can 


be  expressed  (Wilks  (19(')2),  p.  151,  Problem  6. 13]  as 


(L/h) 


(x"  )T(N  -X  )■•  i^  i^ 


(2.2) 


The  p 's  can  be  interpreted  as  the  avc:}\i(jc  rate 

J * ^ 

★ 

over  the  interval  ) : that  is,  p is  the  conditional  prob- 

ability  that  an  item  which  survived  to  time  t^  will  fall  by  time 
t.  . In  terms  of  the  j , we  can  write 

J ► ^ 

X . , “ P,  . + o(h)  . 

J • 1 J * * 

Fquation  (2.2)  provides  us  with  a starting  point  for  a 
Eaves ian  analysis. 


2.1  Prior  Distributions  for  p 


We  shall  confine  the  discussion  in  this  section  to 
the  jth  environmental  condition. 


6 - 


w 
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Snp|>iisi’  that  lh«'  p 


I . I 


I - I . . (I  /h>  . .ir.' 


u indopcndont  beta  random  variables  with  marginal 

densit les 

-I 


■■‘"j.i"  “j.i’ 


* *^*11”^  * i~‘ 


It  is  of  interest  to  note  (cf.  bochner  (197‘>)I  that 
the  above  prior  density  leads  us  to  a -jt'nt'mliai'd  idrid'K! >•: 
dcKsit.  for  the  ^ >\l  . (l./h)  + l ^ • 

Civen  the  N 's  and  the  x ,'s  , the 

J * i J » ^ 


donsityof  . (L/h) 


) is 


a./h) 


a ,+  X , -1 

r ...  J.t  " ...  (n  ) ’ 


* i+N  -X  -1 

(1-p^^^) 


(J.3) 


see  Deuroot  (.1970)  p.  IbO. 


3 . A Bayesian  Kstimation  of  Ordered 
Averap.e  Failure  jbt t_e s 

In  a Bayesian  context.  Condition  (1.1)  leads  us  to  the  require- 
ment that  for  every  fixed  value  of  i , 

P4  4 I.  P«  4 J “ (-3.1) 

J~i«i  J>i 

)Sr  ■ * ^ 

Thus,  our  prior  distributions  on  P^  j Pj  ^ will  have 

to  be  chosen  such  that 

PIp*_j  i 1 p1  1 *’tPj,i  1 P^  P 1 • (3.2'> 

One  way  of  achieving  Condition  (3.2)  is  to  assume  that  the 

* * . - 

parameters  of  the  prior  distributions  of  P^_X  ^ Pj  j satisfy 

- 7 - 
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the  followlnj;  conditions  for  every  fixed  value  of  i : 

«j,i  - “j-1,1 

and  (3.3) 


^j,i  --  J=2.3,...,k. 

For  a proof  of  the  above  statement,  see  Appendix  A. 

In  order  to  be  assured  that  Condition  (3.1)  is  also  satisfied 

with  respect  to  the  posterior  distribiition.s  of  pt  , ^ ‘'nd  p*  , 

j ~1 » 1 J » 1 

it  is  sufficient  (see  Appendix  A)  that,  for  every  fixed  value  of  i 


j.l  J.i  - J-1,1  J-1,1 


and 


(3.4) 


j .1 


~ ’'j,!  - ‘^j-l.i  ”j-l,i 


- X 


,i-l,l  for  j = 2,3, . . . ,k 


In  order  to  make  our  Bayesian  analysis  more  practical,  we  will 
have  to  reduce  tlie  number  of  prior  parameters.  One  way  of  doing  this 
is  to  assume  tliat 


and  that 


M,i 


n for  all  values  of  1 and  j , 


(3.5) 


* 

Thus,  the  prior  distribution  of  each  P,  . > 3 = 1 , 2 , . . . , (I,/h)  , 

J » ^ 

is  a beta  with  parameters  a and  6^  . This  plus  Conditions  (3.3) 
and  (3. A)  lead  us  to  the  following  remarks. 


The  first  part  of  Condition  (3.4)  will  be  satisfied  wiicnever 
’‘i  1 f.  ’‘i  1 “ 2,3,...,k  . That  is,  the  number  of  failures 

in  the  Interval  (t^^t^+h)  under  environment  E|  must  not  hr  tf/'ro/i')' 
i/um  the  number  of  failures  in  ti»e  same  Interval  under  environment 


3 


I - w > 


^ , for  all  valui's  t>f  1 . Horausi*  Ion  (l.'il  will  havi*  to 

bo  aatlsfloii  for  ovorv  fixoit  valiio  o(  J , ( « 1 ,2, . . . , (l./li)  , .1 


roiisounblo  stratoy'.y  Is  to  luivo  N ^ N for  1 » 2,1, 

1*1  I ^ * 


.k  . 


That  Is,  wo  shall  put  moro  Iti'ins  oi\  tost  (lultlallv)  uiulor  t ho  itu'i'o 

• t ho  oav  I rouinont  K j . 

, I “ l.i? (l./h),  tho 


sovoro  ouvlronmoat  f. | j , than  uiulor  tho  onvlronmont  K j . 


Sliu'o  N 


1.1 


^1,1-1  " ’‘l.t-l 


socoiul  pari  ot  I'.oiulltlou  (1,4)  can  bo  wrlttou  as 


1-1, 1+1  - '.1  i.i+i  ' ' 


.k 


( f . ('  ) 


An  Int orprotat ton  of  the  abovo  ooiulltlon  is  that,  for  ovory 
flxod  valuo  ot  1 , tho  nunihor  surviving  at  tho  start  of  tho  t,l-tl)st 
tntorval  plus  tho  I'rlor  paramotor  ^ li>r  t lu-  0110  i |•^>llmc•Ilf  *'j_i  • 


must  bo  smalli'r  than  tho  oorrospoiui  Ing  sum  lor  tho  onvlroniiuMit  f 


Thus,  wholhor  tho  sooond  part  of  tliuuilt  liui  1,1.4)  Is  sat  1st  h\l  or 
not,  ilopoiuls  not  only  on  tho  numbor  of  falluros  In  tho  Ith  Inl.'ival 
aiul  tho  numbor  surviving,  but  also  on  tlu'  valuos  ol  t h«'  prior 
paramotors  bj  and  bj  ^ . Slnoo  tho  numbor  of  tailuros  In  .a 

partloular  Intorval  Is  a funot  Ion  of  tho  sovoritv  ol  tho  onvir.'u- 
mont  and  tho  numbor  on  tost,  and  slnoi‘  has  to  bo  loss  than 

or  oqunl  to  t 


(soo  I’ondltlons  (l.S)  and  t !.(>)),  It  Is  roason- 


ablo  to  havo  b 


I'l 


bj  whonovor  '’i  ''  l.l  • 1‘*.  < !'*• 


valuos  of  tho  prior  paramotors  b]_j  and  b|  aro  indliatlvo 
of  tho  ri'l.'itlvi'  sovi'iltv  ot  tho  onv  I ronmont  a I oonditions 


'.l-l 


and  K , 


Sliuo  K y l',^  ...  > I',  . wo  will  ohooso  tho 


b^'s  suoli  th.it  bj  ' b, 


b|^  , and  tlu-  v.'iluos  ot  tho  bj's 


will  bo  liulloatlvo  of  tho  sovoritv  ot  tho  mw  I ronmont  a 1 conditions 

K,  . 1 - 1...’ k . 

If  tlu'  prior  paramotors  bj  , 1 " l,2,...,l(,  aiul  tho  diita 


, fl./h)  , 


from  tho  accoloratod  life  tost,  x,  , and  N , I " 1,2, 

aro  such  that  I’.omlit  Ion  (1.4)  Is  satisfied  tor  ovorv  flxod  v.iluo  ol 
I , thon  tho  stoch.'isllc  ordorlng  Condition  (,l.l)  will  bo  aut  omat  lea  I 1 v 
sallsflod  with  rospoct  to  tho  posterior  distribution  of 


_ i)  - 


^ 'ft 

(p  ,p  p,  /,/.•>)  • If  the  above  is  not  the  case,  then  we 

will  have  to  pool  the  adjacent  violators  using  the  pooling  procedure 
described  in  the  next  section. 


4.0  The  Pooling  of  Adjacent  Violators 


The  procedure  for  pooling  adjacent  violators  described  here 

is  commonly  used  in  isotonic  vegncasion;  see  Barlow,  Bartholomew, 

Bremner  and  Brunk  (1972).  The  pooling  is  between  the  violators  of 

the  assumed  ordering;  that  is,  whenever  x . . < x , then 

J"-*-,!  J*i 

X and  X . are  pooled. 

j“i» i J » t 

Consider  the  time  interval  l(i-l)h,lh);  by  Condition  (3.4) 


we  require 


V >v  >x 

’'l.l-  2,1-  -’‘j.l  - -*k,l 


“l  * "l.l  - "l.!  - *2  "2,1  - “2,1  i •••  - ®i-l  * "1-1.1  - '‘1-1,1 


*"k.l-'‘k,l 


If  a reversal  occurs,  that  is,  if  either 


’‘j-l.i  ^ "'j.i 


or  if 


^j-1  ^-1,1  " *j-l,l  ""  ^j  *^j,l  " ’'j,!  ’ 

then  we  pool  the  violators  and  replace  them  as  shown  below. 
Replace  both  j j hy 


“^j  ”j-l,i  ^,1  " '‘j-1,1  " ’‘j,!^  • 
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We  now  test  to  see  if  the  new  sequence  Is  properly  ordered,  l.e., 

X,  , > X,  , > Xj,i)  - i(Xj.i_i+  Xj_^)  > 


and 


+ '•’i.i  - -1,1  i •••  i i'^j-i-"  V "j-i.i'"  ”j,r  -j-i.r  -j.i’ 

■ ? * “j-i.i-"  "j.i-  -j-i.i-  -j,i>  i i V\.i-  -k,i 

If  a reversal  exists  in  either  of  the  new  sequences,  then  we 
replace  again  by  appropriate  averages.  Thus,  if 

or  if 

2(3j_i+  6j  + , 

then  we  replace  each  one  of  the  three  by  the  corresponding  average 

3^’‘j-l.i'^  *j,i'^  '‘j+l.i^ 


(in  the  first  sequence)  and 


(in  the  second  sequence). 

We  continue  the  above  procedure  until  all  reversals  in  the 
interval  [(i-l)h,ih)  are  eliminated.  We  use  this  pooling  scheme  for 
each  of  the  (L/h)  time  intervals  to  achieve  the  desired  ordering. 


4.1  Some  General  Comments  Regarding 
the  Pooling  Procedure 

The  following  comments  regarding  the  pooling  of 
adjacent  violators  will  be  helpful. 


I 


1 


- 11  - 


T-375 


1.  An  excessive  amount  of  pooling  occurs  if  th<» 
relationship  specified  by  Condition  (1.1)  is 
incorrect,  or  if  the  environmental  conditions 
Ej,  E2,...,Ej^  are  too  similar  to  each  other. 

2.  If  XjCO'*’)  HO  for  j - l,2,..,k  , then  it 

is  reasonable  to  expect  that  some  pooling  will 
be  necessary  at  the  lower  values  of  i , 
i “ 1,2 (L/h)  . 

For  many  practical  situations,  it  is  reasonable  to 
assume  that  the  failure  distributions  Fj  , j “ l,2,..,k,u  , 
have  increasing  failure  rate  (see  Barlow  and  Proschan  (1975)). 
If  we  wish  to  make  such  an  assumption,  and  Incorporate  it 
into  our  analysis,  then,  for  each  value  of  j , we  must  have 
in  addition  to  Condition  (3.1)  , 


^j.l-l  - ^J.i  ’ 


or  that 


* St  * 


‘■‘•2 


The  above  condition  will  further  complicate  our  pooling 
re;  we  shall  ther« 
increasing  failure  rate. 


procedure;  we  shall  therefore  not  assume  that  the  Fj  are 


5.  The  Posterior  Distributions 
of  the  p*  's 
llii 

It  is  apparent  from  the  discussion  in  Sections  3 and  4 that  the 

* , 

posterior  distributions  of  the  p 's  depend  on  the  outcome  of  our 

j » i 

pooling  procedure.  The  posterior  analysis  is  straightforward  if  no 
pooling  is  necessary,  for  then,  the  posterior  distribution  of 

<pJ.i’p^2---Pj,(L/h)+i> 
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r(a+e.+N  ) 

_J J-li — 


(L/h) 


. a+x  .-1  . 


for  J “ ly2y««»}k.  • 

Under  the  assumption  of  a squared  error  loss  function,  a Bayes 

estimator  of  p , 1 = 1,2, . . • , (L/h)  , j = 1,2,. ..,k  , is  simply 
J > i 

the  posterior  mean;  that  is 


(a+6j+Nj 

(see  DeGroot  (1970)  p.40)  . 

If  pooling  is  necessary,  then  some  or  all  of  the  ^'s  and 

the  (6,+N,  .)'s  will  be  replaced  by  their  appropriately  pooled  averages, 
j J » i 

A ^ 

In  any  case,  the  general  expression  for  p will  be  of  the  form 

j,l 

given  above. 


6.  A Model  for  Extrapolations  to 
Use  Conditions  Stress 

Our  analysis  leads  us  to  an  array  of  Bayes  estimators  of  the 

V 

average  ‘failure  rates  over  intervals  of  length  h , for  each  environ- 
mental icondition.  Because  of  our  pooling  strategy,  the  Bayes  estimators 
will  be.^orrectly  ordered.  That  is,  for  each  fixed  value  of  i , 


'*  ** 

’l,i-P2,i 


P,  , 1 P,  . i •••  1 Pu 


* 

’k,i  ' 


(6.1) 


Given  the  p , s , j = 1,2,.  ..,k  and  i ••  1,2 , . . . , (L/h)  , 

J » ^ 

* * * , 

our  goal  is  to  predict  Py  i » Py  2 ' ” ’^u  (L/h)  ’ average 

failure  rates  over  the  time  intervals  [(i-l)h,ih)  , i “ 1,2, . . . , (L/h)  , 
respectively,  under  the  use  conditions  environment  E . 
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In  the  absence  of  any  knowlodj’,e  about  a relationship  between 
the  average  failure  rates  and  the  values  of  the  various  stresses  whlcli 
constitute  an  environment,  some  form  of  an  assumption  is  essential, 
llils  is  particularly  crucial  if  k , the  number  of  distinct  environ- 
mental conditions  of  interest,  is  small.  If,  however,  k is  large, 
then  a relationship  between  the  average  failure  rates  and  the  stresses 
can  be  empirically  obtained;  this  is  what  is  often  done  in  practice. 

V.lien k is  small,  we  shall  postulate  tlie  following  simple  but 
reasonable  ihi^^ationship  between  the  average  failure  rate  estimates. 

For  some''',  k unkncun  constants  w^  , w^  , . . . , Wj^  ^ , we 
assume  that  for  each  value  of  i , i = 1,2 , . . . , (L/h)  , 


Pk  i • Wq  + Wi  p^  1 i + W2  Pj^  2 i ^ ^ '^k  1 ^1  i 


(6.2) 


The  above  relationship  states  that  the  average  failure  rate 
over  a particular  time  Interval  under  the  environmental  condition 
, is  a weighted  sum  of  the  average  failure  rates  over  the  same 
time  Interval  under  the  conditions  2 ’***’  ^1  ‘ This  is 

reminiscent  of  an  autoregressive  process  of  order  k-1  which  has  found 
useful  applications  in  forecasting  (see  Box  and  Jenkins  (1975)). 

In  order  to  make  the  above  relationship  more  meaningful  we 
shall  require  that  the  environmental  conditions  Ej^<  ^ K-) 

increase  in  magnitude  of  severity  by  the  same  fixed  amount.  For 
example,  if  represents  a single  stress,  say  a voltage  stress  , 

then  we  shall  require  that  Vj  = C , where  C is  some  suitable 

constant,  j = l,2,***,k  . 

Since  Equation  (6.2)  holds  for  i = 1,2 , • • • , (L/h)  , the  least 
squares  estimators  of  Wq  , w^^  ’’'*’'^k-l  obtained  in  a routine 

manner.  These  estimators  are  denoted  by  , • • • , • 
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If  the  Increase  in  the  magnitude  of  the  severity  of  the 

environment  from  E to  E,  is  the  same  as  that  from  E,  to  E . , , 

u nt  j J-1 

J ■ 2,3,**’,k  , then  an  estimator  of  the  average  failure  rate  under 

E as  given  by 
u 


Pu.i  ■ '^0  Vk.l+ 


Vl.i  P2.i 


(6.3) 


(L/h) 


for  i - l,2...-,(L/h)  . and  (L/h)+i  ^ P^.i  * 

If  the  increase  in  magnitude  of  the  severity  of  the  environment 
from  E to  E,  is  two  times  that  from  E to  E ^ , j = 2 , 3 , . . . ,k  , 

UK  J j~l 

then  we  iterate  upon  Equation  (6.3)  one  more  time  to  obtain  the  desired 
result.  Thus,  in  principle,  we  can  Iterate  upon  Equation  (6.3)  as  many 
times  as  is  necessary,  depending  upon  the  separation  between  E and 

\ • 

* 

By  the  definition  of  p . , the  probability  of  an  item  surviving 
to  time  t (assumed  a multiple  of  h ) under  environment  E^  is 

_ * * 

F (t  ) = n (1-p.  ,)  . 


i=l 


"u,l' 


Thus,  we  have  the  following  as  an  estimator  of  ^^(6  ) •’ 


1 * tVh 

Fu(t  ) = (i-p^^^)  , 


where  the  p . are  given  by  (6.3). 

U » 1 

The  properties  of  this  estimator  have  not  been  studied.  The 
estimation  procedure  has  been  proposed  as  a practical  method  for 
providing  answers  in  the  difficult  area  of  accelerated  testing. 
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Appendix  A 


In  this  appendix  we  determine  sufficient  conditions  under  which 
two  beta  random  variables  are  stochastically  ordered.  The  results  may 
already  be  known;  however,  they  are  included  for  completeness. 


The  notation  and  terminology  are  taken  from  Barlow  and  Proschan 

(1975). 


We  first  note  that  a gamma  density  f(x,a) 


a-1  “X 
X e 

r(a) 


is  TP^ 


a-1 


in  (x,a)  . Therefore,  F(x,a)  = / / y e ^dx  is  TP^  in  (x,a)  since 


F(x,a)  “ J f(y,a)  K(x,y)dy,  where  H(x,y)  = 1 for  y ^ x and 


H(x,y)  » 0 for  y < x . (Note  that 


H(Xj^,yj^)  H(Xj^,y2) 

n(x2,yj^)  H(x2,y2) 


> 0 for 


*1  ^ *2  ’ ^1  ^ ^2  ' H(x,y)  is  TP2  , and  so  the  composition  of 

f,  H is  TP,  .) 


Since  F(x,a)  is  TP2  , then  for  x^^  < X2  , < 02  , 


F(Xj^,a^) 

P(xj^,a2) 

1 0 . 

Choose  X,  » 0 . Then 

1 

1 

> 0 

F(x2,a^) 

F(x2,a2) 

1 

F(x,oij^) 

F(x,a2) 

for  0 < X 

< * , 

<02  . 

Thus,  F(x,a2)  ^ F(x,aj^) 

for 

<a2 

and  X > 0 . 
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We  have  shown: 


1.  Lemma.  Let  be  a Ramma  r.v.  with  shape  parameter 

a . Then  X is  increasing  stochastically  in  a . 


2.  Proposition.  Let  Y 


have  density 


1 (,~y  1 = 12 


and  Y , Y be  independent.  Then  X = v — Z'v — has  beta 

^ “2  “l*“2  VS 


r(a^+a„)  ®i“^  °‘2~^ 

density  f(a^.a,,x)  = x (1_^)  0 <1  . 

(Hogg  & Craig,  1970,  p.  13A). 


3.  Theorem.  Let  X have  beta  density  f(a  ,a^,x)  . 

“r«2  ^ “ 

Then  X is  increasing  stochastically  in  a and  decreasing 

ai.«2  1 

stochastically  in  a2  . 

Y Y 

“l  "'l 

Proof.  Write  X * — —v — » ^ t = 1. — ~r  v~~  > 

a, .a-  Y + Y n, »nA  V + Y , ’ 

12  “i  *^2  ^ ^ “l 

St 

where  a.,  < a'  . For  fixed  Y ,X  >X  , _ j 

2 2 a^^  ^1’‘‘'2  ~ '^1’°'2  ' unconditioning 


'■^’“2  " "^a^.a-  «2  ""2  • 

St 

By  similar  reasoning,  we  can  show  X < » for 


“l’”2  ~ “l’"2 


Of  < «;  . 
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